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Abstract 

We study the 0(AQCD/'"^fe) corrections to the singly and doubly differential hadronic invariant 
mass spectra dV/dsn and dV / dsndq^ \n B ^ X^ii^e decays, and discuss the implications for 
the extraction of the CKM matrix element \Vub\- Using simple models for the subleading shape 
functions, the effects of subleading operators are estimated to be at the few percent level for 
experimentally relevant cuts. The subleading corrections proportional to the leading shape function 
are larger, but largely cancel in the relation between the hadronic invariant mass spectrum and 
the photon spectrum in B ^ Xg^. We also discuss the applicability of the usual prescription 
of convoluting the partonic level rate with the leading light-cone wavefunction of the b quark to 
subleading order. 
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I. INTRODUCTION 



The CKM parameter \Vub\ is of phenomenological interest both because it is a basic 
parameter of the Standard Model and because of the role it plays in precision studies of 
CP violation in the B meson system. Currently, the theoretically cleanest determinations 
of \Vub\ come from inclusive semileptonic decays, which are not sensitive to the details of 
hadronization. 

For sufficiently inclusive observables, inclusive decay rates may be written as an expansion 
in local operators Q|. The leading order result corresponds to the decay of a free b quark 
to quarks and gluons, while the subleading corrections, proportional to powers of AQCD/''^fe, 
describe the deviations from the parton model. Up to ©(Aq^^p/m^), only two operators 
arise, 

Ai = J—{B\K{tD)^K\B), A2(/i) = -^{B\Ka^''G^,K\B). (1) 

The B — B* mass splitting determines A2(7Ti;,) ~ 0.12GeV^, while a recent fit to moments 
of the charged lepton spectrum in semileptonic b ^ c decay obtained 

m^^ = 4.82 ± 0.07s iO.llTGeV, Ai = -0.25 ± 0.025^ ± O.OS^y ± 0.14t GeV^ (2) 

where ml^ is the short- distance "IS mass" of the b quark . (Moments of other spectra 
give similar results 0,0].) These uncertainties correspond to an uncertainty of ~ 5% in the 
relation between \Vub\ and the inclusive B Xuii'e width 0,0|- 

Unfortunately, the semileptonic b ^ u decay rate is difficult to measure experimentally, 
because of the large background from charmed final states. As a result, there has been much 
theoretical and experimental interest in the decay rate in restricted regions of phase space 
where the charm background is absent. Of particular interest have been the large lepton 
energy region, Ei > {m'^ — m'j^)/2mB, the low hadronic invariant mass region, mx = < 

TTlD 



(jthe large lepton invariant mass region > (m^ — 0], and combinations of 

these 10|. The charged lepton cut is the easiest to implement experimentally, while the 



ladronic mass cut has the advantage that it contains roughly 80% of the semileptonic rate 



IjJ. However, in both cases the kinematic cuts constrain the final hadronic state to consist 
of energetic, low-invariant mass hadrons, and the local OPE breaks down (this is not the 
case for the large region or for appropriately chosen mixed cuts). In this case, the relevant 
spectrum is determined at leading order in AqcB/^b by the light-cone distribution function 
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of the b quark in the meson 12, ll3 1 



{B\b6{u + in- D) h\B) 
2mB 



(3) 



where is a hght-hke vector, and hatted variables are normahzed to m^: Z)^ = D^/rrif,} 
f{uj) is often referred to as the shape function, and corresponds to resumming an infinite se- 
ries of local operators in the usual OPE. The physical spectra are determined by convoluting 
the shape function with the appropriate kinematic functions: 



1 ir 



To dEi 
1 dT 



4 / e{l-2Ei-u)f{uj) duj + 
2s|^(3cj - 2sh) 



SH)f{^ — A) duj + . . . 



(4) 
(5) 



To dsH 

where 1 - 2Ei < Aqqd /rrib, Sh ^ Aqcd /rrib and A = niB - mi,. 

Since f{uj) also determines the shape of the photon spectrum in B ^ Xg'j at leading 
order, 



1 dT 



{B ^ X,7) = 2/(1 - 2E,) + 



(6) 



r§ dE^ 

there has been much interest in extracting f{uj) from radiative B decay and applying it to 
semileptonic decay. However, the relations (jlHSl) hold only at tree level and at leading order 
in Aqcd/'"^;), so a precision determination of \Vub\ requires an understanding of the size of 
the corrections. Radiative corrections were considered in ll^ 13 . Il^ lisl . while ©(AqcD/'^ft) 
corrections have been studied more recently in 1^, ^ lisi l^j. In |l(|. the nonlocal dis- 
tribution functions arising at subleading order were enumerated, and their contribution to 
B decay was studied. In l3|, the corresponding corrections to the lepton endpoint 

spectrum in i? — > X^fz/^ decay were studied, and it was shown that these effects were po- 
tentially large. Similar results were obtained in where the sub-subleading contribution 
from annihilation graphs was also shown to be large. In this paper, we study the sublead- 
ing corrections to the hadronic invariant mass spectrum in semileptonic 6 — > m decay, and 
estimate the theoretical uncertainties introduced by these terms. In addition, we present 
results for the doubly differential spectrum dV /dsudq^ at leading and subleading order. 



^ Because in our definition of /(w) its argument is dimensionless, /(w) differs by a factor of rrib from the 
usual definitions in tlie literature. 
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II. MATCHING CALCULATION 



A. The full theory spectrum 

In the shape function region the final hadronic state has large energy but small invariant 
mass, and so its momentum lies close to the light-cone. It is therefore convenient to introduce 
two light-like vectors and related to the velocity of the heavy meson by = 
]^{n^ + n^), and satisfying 

= = 0, V ■ n = V ■ n = 1, n ■ n = 2. (7) 

In the frame in which the B meson is at rest, these vectors are given by = (1,0,0, 1), 
= (1,0,0,-1) and v'^ = (1,0,0,0). The projection of an arbitrary four-vector a" onto 
the directions which are perpendicular to the light-cone is given by a" = gx^ap, where 

Choosing our axes such that the momentum transfer to the leptons q is in the — n direction, 
we can write = ■ q n'^ + ■ q n^, the decay rate takes a particularly simple form in 
terms of the variables n ■ q and n ■ q: 

dT{B ^ XJue) = 967r Tq W^,^Lf"'{n -q-n- qfe{n ■ q)e{n -q-n- q)dn ■ q dn ■ q (9) 

where 

° 1927r3 ^ ' 

The hadron tensor W^" is defined by 

,y.^^_ito('_,:|Ae---<M|Ma), (11) 

where the weak current is Jj^ = ^7^(1 — '~f^)b, while the lepton tensor is 

L^^ ^ J dU2iq;pe,p.)TTy>,7%YPL] = - (12) 

andPL = i(l-75). 

To calculate the hadronic invariant mass spectrum we switch to the variables {suiq^)- 
These are related to the variables in Eq. Q by 

Sh = {rrtB - n ■ q){mB - n ■ q) 

= {nib + A — n ■ q){mb + A — n ■ q) (13) 
= n ■ q n ■ q (14) 



and 

= V((™. + A)^ + ,-.„)-4(,„. + A)V^. (15) 

Here A = rriB — rrih is the difference between the B meson mass and the b quark mass. It is 
(9(Aqcd) and has an expansion in terms of HQET parameters 

A=A-^if^ + ..-. (16) 

Since A simply enters in the definition of sh, it is unrelated to the 1/m expansion in the 
OPE, so we will not expand it via Eq. (fTI)|) . With this change of variables, we define the 
correlator T(s//, g^) by 

^ ^ ■{msH,q')\B). (17) 



To dsH dq^ 2mB 

In Ref. [l6| a nonlocal expansion was performed for the hadron tensor W^'^, based on the 
power counting 

{rribV — q) ■ n = nib — q ■ n 0{nib), 

{nibV -q)-n = nib-q-nr^ C(Aqcd), (18) 
k'' ~ O(Aqcd). 

where the heavy quark momentum is defined as Pb — ^bV^ + k^- However, the limits of 
phase space integration in Eq. include regions of phase space where this power counting is 
violated. Hence, to keep our power counting consistent, we do not perform a nonlocal OPE 
for W^" , but rather for T{sh, q^)- In these variables, the shape function region corresponds 
to the region of low invariant mass, 

Sh ~ C(AQCDmb). (19) 

Since A ~ Aqcd and k^^ ~ Aqcd, expanding the light quark propagator in powers of 
AqcD/^^b gives at leading order 

pI 2nib[sH-{k-n-k){l-(p)] 

(where A = A/nib). Since both terms in the denominator are 0(Aqcd/''^6), T{sh, q^) cannot 
be expanded in powers of k^ and matched onto local operators (unless we also are restricted 
to large g^, such that 1 — -C 1, in which case the second term in the denominator 
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is subleading, and a local OPE may be performed [9|, llQ|). Instead, the OPE takes the 
schematic form 

T{sh, q^) = Y. [ sh, f)Oniuj)dLU (21) 

where the (9„(c(j)'s are bilocal operators in which the two points are separated along the 
light cone. 



B. Nonlocal operators 



In Refs. 



[7 



it was shown that up to subleading order in Aqco/'^^b, the following 



operators were required in the OPE 

Oq{uj) = h^S{uj + in ■ D)hy 

h,^{iD^',5{uj + in ■ D)}K 
hy[iD^,6{uj + in- D)]hy 



(22) 



dujiduj2 6{iUi, UJ2; uj)hy5{in ■ D+uj2)g^ {iDf^, iD^}5{in ■ D+uJi)h 



Oi{uj) = — J duidu2 5{ui,U2;uj)hyS{in ■ D+U2)e±'[iDfj,,iD,y]S{in ■ D+uJi)hj; , 

Pq{uj) = hy5{LU + in ■ D)'^'^'^^hy 
pf (cu) = Kiib", 5{uj + in ■ b)}-i'^^^K 
PPiuj) = KliD^", 6{uj + in ■ D)]-f'^^5K 

P^{uj) = J duoiduj2 5{uoi, ijj2\ uj)hy5{in ■ D+uj2)g'^{iD^, iDy}5{in ■ _D +co'i) 7*^75 
P4 (cu) = — j dujiduj2 6{uJi, 102] uj)h.u5{in ■ D+uj2)e'^[iD^,iD^]5{in ■ D+uJi)'j'^'j^hi 

where the h^^s are heavy quark fields in HQET, and we have defined 

6{uj — ui) — 6{lj — UJ2) 



6{uji,uj2;u}) 



and 



0/3 



Ui - LO2 



(23) 



(24) 



These definitions differ slightly from the definitions in Refs. [16|, ll7[ , because we have chosen 
to normalize all momenta to m?,, to keep the resulting formulas simpler. 

It is convenient to calculate the matching conditions onto a slightly different set of oper- 
ators, defined in terms of full QCD h quark fields: 



Qo('^,r) = U{in-T)\uj)Vh 



(25) 



<5o(a;,r) 




5{u) + n ■ k)T 




2k"5{uj + n ■ k)T 





Q^(<^,r) 





a, a 

grV5_(n-^)r 



Q3('-',r) 



a, a 



<34('^>r) 





\ n - e 




^-("■(^1+^2)) 
n • (ii + 4) 



FIG. 1: Feynman rules for the operators Qi{uj,T) in n • ^ = gauge. We have defined S±{x) 
6{ijj + n ■ k + x) ± 6{uj + n ■ k) . 



Q3(w,r) 
Qi{uj,T) 

We have defined 



b[iV'',6{in-V + uj)}Tb, 



i'D^',5{in-'D + uj] 



r6, 



duJiduj2 S{uJi, UJ2] uj)b6{in ■ V+uj2)g'^{i'D^, iVy]5{in ■ V+uijVh, 
dujiduj2 S{uJi, UJ2; u!)b6{in ■ 'D+uj2)e'^[iV^, iVy\5{in ■ V+uJi)Th . 



(26) 



so that iV^ acting on the h fields just bring down factors of the residual momentum k^. The 
Feynman rules for the Oj's and Pj's are given in The rules for the Qi's are given in 

n ■ A = Q gauge in Fig. [H where we have defined 



5±{x) = S{u! + n ■ k + x) ± 6{uj + n ■ k). 



(27) 



It is simpler to match onto the Qj's initially since this matching does not require us to relate 
the QCD quark fields to HQET quark fields. However, because the additional symmetries 



of HQET reduce the number of independent functions needed to parametrize the matrix 
elements, it is convenient to then express the QiS in terms of the Oj's and Pj's. For an 
arbitrary Dirac structure F we have 

b5(u} + in ■ V)rb ^h^(l + ^ ) S(cu + in ■ D)r ( 1 + ^ ) + ■ ■ ■ (28) 

V 2m6y \ 2mbJ 

= ^IV[FP+]Oo(a;) - ilV[F5,]P;(c^) 
+ ^ (Tr[7,FP+] + IV[F7aP+]) O^iu) 
+ ^ (I¥[7,FP+] - Tr[F7,P+]) 0^{u;) 

- ^(IY[7,F5,]+Tr[F7,5,])A'''(a;) 

- ^ (IV[7,Fs,] - IV[F7,s,]) P.^^'ico) + ■■■ 

where 

P+ = i(l + ^) and = P+7''75P+ (29) 

and we have used the fact that 

h^rh^ = ^Tr[rP+]h^h^ - ^Tr[Fs^] ^^7^75/1^ . (30) 

For our purposes, we will only need the case F = jaPL, which allows us to write 

Qo{uj, -faPl) = ^VaOo{uj) - ^(fir^^ - VaVr,)PSiuj) (31) 

+ ^i(^aX,^pV^P2'^{i^) H 

where the first line gives the leading order relation and subsequent lines contain the sub- 
leading correction. 

Similar relations may be derived for the subleading operators, though in these cases it is 
not necessary to consider the subleading terms in the relation between the QCD operator 
and the HQET operator, such terms being of higher order overall. Thus we have 

Ql{oj, 7,P^) = Iv^O^ioo) - \{g.,„ - v^v,)P^\uj) + ■■■ 

Q^iuj, j^Pl) = lv,0^{uj) - \{g^„ - v,v,)Pr{'^) + ■■■ (32) 

Qsii^naPz) = lv^03{u;) - ^{gr,a - Vr,v^)PS{u;) H 

QA{<^,laPL) = \VaOA{uj) - \{gna - VrjVa)P2{(^) H ■ 



The leadin g an d subleading operators can then be completely parametrized in terms of 
five functions Il6ll: 



(fi|Oo(^) 


\B) 




\B) 




\B) 




\B) 


{B\P2{uj) 


\B) 



2mB f{uj) + 



t{uj)- 



(33) 



(once again, unlike in these are defined here in terms of dimensionless arguments). The 
matrix elements of the other operators vanish. 



C. Matching Conditions 

The Wilson coefficients Ci{u!) of the operators in (PT|) are obtained by taking partonic 
matrix elements of both sides of the OPE. In particular we take zero-, one-, and two- 
gluon matrix elements, which corresponds to calculating the imaginary parts of the full- 
theory forward-scattering diagrams in Figure El multiplying by the lepton tensor L'^^ and 
appropriate phase space factors and matching them onto linear combinations of the effective 
diagrams. (The matching conditions may be completely determined from just the zero-gluon 
and one-gluon matrix elements, but we have calculated the rest as a check of the results.) 

The lepton tensor has the expansion 



[nf'n'' + q^in^n" + n^'rf) + q^n^'rf - Af g^""} (34) 
2^ [n^'n'' - q^n^'n"} + ■ ■ ■ . 



4871 



247r (1 - g2 

(where we have used the decomposition q^ = n ■ + n ■ qn^ j1\ while the phase space 

factors give 

in - q — n- g)^ 



((l + A)2 + g2_s^)2_4(i+A)V 

= (i-mmi-a 



{n ■ q — n ■ q)6{n ■ q) 



1 — 



a, a 



(a) 



(b) 




cti, a ct2, b 




(c) 



FIG. 2: Full-theory forward scattering diagrams. 



The zero-gluon diagram in Figure a) gives the amphtude 

Pi 

Taking the imaginary part of this amphtude gives 



(36) 



TT 

= ^rfiYPL [(1 - f)Sihin ■ k)) 



+ 



sh 



(37) 



(A(l - q^) - snr) - (1 - r)ki 6'{h{n ■ k)) 



+Yhl"PLKh{n ■ k)) + 



where we have expanded the amphtude to subleading order using (fT^ and we have simphfied 
the expression by integrating by parts. The function h{x) appearing in ()37|) is 



h{x) = Sh -{A- x)(l - 



(38) 



Multiplying this result by the lepton tensor and phase space factors and expanding 
to subleading order we find 



{b\TisH,f)\b) = Y. duc:{u;,SH,e){b\Qn{oo,rPLm 



(39) 
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where 



Co^(^, SH, (f) = 4(1 - g2)^(2g V + n'^)d{e)d{l - e)6{h{-u;)) 
+8(1 - g')((A(l - 3f) - ujf{3 - g2)K 

-(A + uj{2 - f))n'')9{qy{l - f)Sih{-co)) 
+4An''6{f)6{h{-uj)) + --- (40) 
CT (^, Sh, f) = 4g^(l - f)gr0{f)9{l - f)S{h{-uj)) + ■■■ 
C^icu, SH, (f) = -2(1 - g')(2gV + n'')e{e)e{l - f)5{hi-u;)) + ■■■. 

In order to determine the other matching coefficients, we calculate the one-gluon ampli- 
tude in Figure |2fb). Defining £ to be the incoming gluon momentum, we have 

where [a, a) are, respectively, the Lorentz and colour indices of the gluon field. 

Taking into account the two cuts which result from taking Im[y4i] and scaling the gluon 
momentum as ~ C(Aqcd), we obtain, after expanding to leading order in n-A = gauge, 

- ^lm[A^] = -^-^r i2ll64n ■ i) + 2/i(2fc + £)° (^-^^] (42) 



_(n ■ 
n-i 



+ 2^ef 7/3755- (n ■ ^) + 2^ef £\^^75 f^^^l I YPl + 



where, in analogy with (P7j). we have defined 

S±{x) = 6{h{n ■k + x))± 6{h{n ■ k)). (43) 

Again, multiplying by the lepton tensor and phase space factors gives 

{h\T{sH,e)\hg)=Y. I duj C:iu;jH,e){b\Qniu;,YPL)\b9). (44) 

Part of ()42|) is reproduced by combining the Wilson coefficients (pn|) determined earlier with 
the one-gluon Feynman rules for (5i_3(co', 7'^Pl), while the remainder corresponds to matrix 
elements of Q2,Aiy>,l'^ Pl) with the coefficients 

CT (^, SH, f) = 4g2(l - e)zeTe{qy{l - q')6{h{-u;)) + ■■■ 
Cl{uj, SH, f) = -2(1 - g^)(2gV + n'')9{f)e{l - f)6{h{-u;)) + ■■■ . (45) 
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The final matrix element to evaluate is the two-gluon amplitude, Fig.|21^c). The amplitude 

is 

so that after cutting the diagrams and expanding to leading order, again in n • A = gauge, 
we obtain 

The two gluon matrix element of T{sh, (f) agrees with the results of and ()45|1 for C3 
and C4; hence, no new operators are required, as expected. 

Integrating these expressions over we obtain the OPE for dT/dsn 
1 dT 



T^dsH 2mB — J-oo 



1 r°° 

— Y^ dujC:{uj,SH){B\Qn{u,^„PL)\B) (48) 



where 



C^{uj, sh) = \, ^-^0 (sh) 9(A + uj-sh 

[A + ujp ^ 

+ .^^"^ {{s^uj + sh{uj^ + 4a;A + sA^) - 2(A + ujf)n'' 
(A + cj)^ 

-(A + uj){A^ + 2uA + uj{sh + oj))n''^ 9 (sh) 9 (^A + uj ~ sh] 
+ 4An'"5 (A + - Sh) + ■ ■ ■ (49) 

Cf [uj, Sh) = -7 9 [sh) 9 [A + u - Sh) ^ 

[A + uj)'^ ^ ' 

{uj, Sh) = -J 9 {sh) 9 [A + uj - Sh) ^ 

(A + cj)3 ^ ^ 

- ^ 23^(2(3^ - A - uj)n'^ - (A + ^)n-) ^ ^ . . x 

C3 (cj, Sh) = -J 9 [Sh) 9[A + uj - Sh) ^ 

(A + uj)-^ ^ ^ 

^ . 2sh(2(sh- A-cj)r2<^- (A + cu)n'^)^_ ,^/;^ N 

Cl{uj, Sh) = ' \ ' '-^9 (sh) 9(A + uj-sh)+---. 

[A + uj)-^ ^ ' 



Finally, relating the Qj's to the Oj's and Pj's via (j3Tj) and (j32|l and taking the matrix 
elements (jSBI)? we obtain the expression for the hadronic invariant mass spectrum: 

1 rfr r , f 2g|,(3^-2gH)/(^-A) ^^. . 

- ' dt^ <j — 9 [Sh) 9[uj- Sh) (50) 



Fq (iSH J-oo 1 UJ^ 
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+ 2Af{u;-A)9{sH)S{u;-SH) 

2sh(4s^(c<; - A) + SHUji7A - cu) - 6uj^)fiuj - A) 



+ 



+ 



+ 



s%{3uj - 2sH)t{u - A) 2s//(3cj - 2sh)G2{uj - A) 



2sh{2s% +UJSH - 2uj^)hi{uj - A) 



2sh{2sh-uj)H2{uo- A) 



e(s„)e(u-s 



H 



+ 



Eq. (jHnil is the principal result of this paper. It may be checked for consistency with the 
result obtained via the local OPE by ex pan ding the matrix elements of the operators (j22p 
such that in ■ D (9(Aqcd)- This gives jl6| 

(cj) H 



/M 






- 5'(- •) \ ^^5 
Sm^" ' 6mf 


hi{uj) 


— 5'{uj) + 5" 
ml 2ml 




2Ai ^, , , 
--,mf^^^^--- 




= + . . . 

mi 


tiu) 


Ai + 3A2 . 
mi 



ISmf 



(51) 



r 



2ml 



where each term in the expansion is of the same order in the shape function region, but 
the terms indicated by ellipses are higher order in the local OPE. The Ai 2 parameters are 
defined in ([T)) and the pi^2 parameters are defined by 



ZmB o 

-^{B\hJDJD^iDpS5K\B) = ^ie^o.psV' v ^p2. 
ZmB Z 



(52) 



When substituted into the spectrum and integrated over u we obtain to subleading 
order 

1 rfPiocal i2sl{3A~2sH) 



To dsH 



-e{A - SH)0irsH) 
+A, i'JM^^i^^oiA - ,M + ^A^. - A)^ 



(53) 



3A6 



3A2 
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+ ^^"^("/;-^) g(A - 5^)0(5^,) + 2A6{SH - A) 

+ Ai ^ ^-^ A - SH)9isH) + -^S{sh - A 

V 3A5 3A 7 

+ A2 i— ^^(A - SH)e{sH) - jS{sH - A) j 

„ /4sh(20s2^-33As^^ + 3A2)^ . ^ ^^^^ ^ 6 A 
+ Pi ( ^^(A - SH)eisH) + ^5{sH - A)J 

. /4si^(10s|^ + 3Asi^-3A2)^ . „ ^^^^ ^ 8 ;>A 
+ P2 ( ^ ^^-^^ ^^(A - sh)9{sh) - -^6{sH - A)J + ■ ■ • 

where the terms in curly brackets are the leading order result, and the other terms are the 
subleading order correction. 

The local OPE spectrum can be obtained from the double-differential spectrum 
dr/dsodEo presented in and Q- 

After changing variables to {sh,Eq) and expand- 
ing in powers of AqcB/mb (treating sh as order Aqco^fe), performing the Eq integral we 
obtain the local OPE for dT/dsn, which exactly reproduces the result 



III. RELATION TO PREVIOUS WORK 

At leading order in l/rrib, the effects of the distribution function f{uj) may be simply 
included by replacing in the tree-level partonic rate 

nib ml = mb(l — a;) (54) 

and then convoluting the differential rate dV with the distribution function f{uj) jl^ . 

dV = j dV^'''^°^\m,-.mlf{uj) du. (55) 

Because of the leading factor of ml in the rate (|Tn|) . this prescription leads to large subleading 
corrections if the factor of ml is included in the replacement ()54|) . 

In Ref. this prescription was applied to the sh spectrum, although the mf term was 
not included in the replacement. This is perfectly consistent at leading order, but since other 
subleading effects were introduced in Ref. jll| by the replacement (|5l| . it is instructive to 
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compare our result (fSUj) with the results of Ref. [ill, expanded consistently to subleading 
order in l/rub. At leading order, the results are identical:^ 

1 dr(o) 



/ A{sH,uj)f{uj - A) du (56) 

^0 



To dsH Jo 
where 

A{sh, uj) = ^-9{lu - sh) (57) 

UJ 

At subleading order, the relevant terms in Eq. (j5n|l may be written as 

1 dr^^^ r°° 

— -p-= 6A{sH,uj)f{uj-A)duj + ... (58) 
1 dsH Jo 

where the ellipses denote subleading shape functions, the effects of which cannot be re- 
produced by the prescription ()55j) . We will refer to these corrections as true subleading 
corrections, and the terms arising from 6A[sh,^^) as kinematic correction. The function 
6A{u, Sh) is 

6A{lu, Sh) = ^ ^ -Oiuj - Sh) + 2A5{uj - sh) 

2sH{Ssjj{A -u) + 3shuj{5u - 3A) - 6u') ^ . , 

= 7 \- 2Ad{uj - Sh) 

lQsl{2sH-3uj){uj~A) 2sl(2sH ~ uj){lo - A) 
+ + 



The second line of Eq. fl3^ agrees with the expansion of the results of Ref. [ll| to subleading 
order. The first term in the third line agrees with the expansion if the mf factor is also 
included in the convolution. Finally, the last term in Eq. ()59p arises from the expansion of the 
quark fields in terms of HQET fields in the relation (f^Hj) . Thus, we see that to be consistent 
to subleading order, one must include the term in the replacement ((SH). However, like 
the subleading shape functions, the subleading effects arising from the expansion of the 
quark fields cannot be reproduced by this procedure. 



The relative sizes of each of the terms in 
one-parameter model for f{uj) introduced in |l 



is plotted in Fig. |2l using the simple 



/^od(^) = 4t(1 + ^/A)'e-^(i+-/^)'^(l + uj/A) (60) 



^ In Ref. the upper limit of integration is a; = y/sn; however, the difference is higher order. In addition, 
the region uj ^ ^fsu 3> sh is outside the region of support of the shape function, and so is expected to 
be suppressed. 
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FIG. 3: Plot of the kinematic corrections to the hadronic invariant mass spectrum, Eq. (|58j) . The 
dashed Une is the leading order result (|56() . while the solid line includes the full set of kinetic 
corrections. The dotted line corresponds to the expansion of the results of 11] to subleading order, 
while the dot-dashed line also includes the contribution from the term. The difference between 
the dot-dashed and solid curves is due to the expansion of the heavy quark spinors. 

and with A = 0.1. Numerically, the most important of these corrections corresponds to 
smearing the ml term, while the correction from expanding the quark fields is quite small. 

However, such large corrections may be misleading, since if they are universal they may 
simply be absorbed in a redefinition of the leading order shape function. Instead, one should 
look at the corresponding relation between the hadronic invariant mass spectrum and the 
B — s> Xs7 photon energy spectrum. One might expect that the effect of convoluting the 
ml term would cancel in the relation, since both rates are proportional to m^. However, in 
the B Xs'~f spectrum only three powers of m^ come from the kinematics, while two arise 
from the factor of in the Wilson coefficient of O7, and hence for this rate one should only 
convolute three powers of mi,. This may be verified by writing the results of Ref. |3| as 
1 dT 



Tq dx 



{B X,7) = /(I - x) - 3(1 - x)/(l - x) + (1 - x)/(l - x) + . . . 



where once again the dots denote additional form factors, and the partonic rate is 



J- 



G'p\VtbVt:\'a\Cf\'ml 
327r4 



(61) 



(62) 



In the expression (|UT|). the second term corresponds to smearing three powers of in the 
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rate, while the third third term arises from the expansion of the quark fields. Thus, there is 
an incomplete cancellation of the kinematic corrections between the two spectra. 



IV. PHENOMENOLOGY 



A. The B Xu^u hadronic invariant mass spectrum and the B Xsj photon 
energy spectrum 

As discussed in the previous section, there are large kinematic corrections to the leading 
order results, largely due to the ml term in the rate. However, these are reduced in the 
relation between the hadronic invariant mass spectrum and the B Xs'~f photon energy 
spectrum. Similarly, the T-product t{x) is universal for all processes involving B meson 
decays (it only differentiates between B and B*, D and D* decays) and so its effects similarly 
cancel. Hence, it is useful to express the hadronic invariant mass spectrum in terms of the 
experimentally measurable B Xg'-f photon energy spectrum. 

?he B photon energy spectrum is given at tree level to subleading order in 1/ rrih 



by 



3 



1 dr 



2F(1 - 2E^ 



(63) 



where 



Fix) = fix) + 



hUx) + 



tix) 



2xfix) -G2ix) + H2ix) 



+ 



(64) 



(Note that at tree level only the operator O7 contributes. At one loop, effects of other 
operators must be included 

1 dV 



2l\). Substituting this into Eq. (jKn|l gives 
rftu {(A(u;, sh) + 6Aiuj, SH))Fiuj - A) + 5F(u;, sh, A)} du 



Tq dSH JO 

where AiujSn) and 6AiLJ,SH) are defined in and (j^Hl), and 

2sHi2sH — 3c<j)(si^ — u 



(65) 



6Ficu,SH,A) 



-G2(cu- A) 



^Asni2sn + uj)isn-u)^^^^_^^ 



, 2sui2s\-^usn^^'') ^ , 
H H2{uj - A) 



diuj - Sh) 



(66) 
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contains the subleading shape functions. (Note that the dependence on the T-product t{x 
drops out of this relation.) 

To extract |Kib|, we are interested in the integrated rate 

rn dV 



up to a maximum value s^. The integrated rate for B X^iu^ is free of backgrounds from 
B Xcii^i for s'fj < mj^, although because of experimental resolution the experimental 
cut is typically somewhat lower: a recent BABAR measurement j^l used sjj = (1.55 GeV)^. 
From Eq. (jHS)), we have 



(68) 



where 



6A{uj, s^) 



e{sH -uj) + e{uj 



8(A -cj) 



UJ) + 



2(s^)2(s^(A + 5cj) - 9cj2 



[UJ 



SF{lj, s'h) = -- {G2{uj - A) + 2hi{uj - A) + H2{lo - A) 
(s^)2(3(S^)2-10cJs^ + 9cJ2^ 



UJ 



+ 



3uj^ 



+ 



2(g^)'(3(g^)'-2cus^-3cu' 



3cj^ 



-G2{uj-A) 



-hi{uj-A) 
■H2{uj - A) U[u 



(69) 



Note that the upper limit of integration in u corresponds to a photon energy = l+A—u < 
0; however, as discussed earlier, this region is outside the region of support of the shape 
function, and its contribution should be highly suppressed. Thus, in the relation between 
the spectra we may set the lower limit on to zero. 

Comparing the two forms for the integrated spectrum ^sh{^h) f!63|) and (j68p . we can 
isolate the CKM parameter \Vub\'- 



IK 



ub 



\VtbVtl\ 



1 - -6T 



6a\C- 



■off 12 \ 1/2 



7 



TT 



X 



dT 



ds 



dSH 



H 



. . . ^ dV 

{A{mB - 2E^, Sh) + 6A{mB - 2E^, SH))^-dE^ 

dE^ 



1/2 



(70) 
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where we have defined 



6F{uj, s'jj) du 

STs, (5^) = : (71) 

A{uj,SH)f{uj - A) duo 



which contains the effects of the new subleading distribution functions. For comparison 
purposes, we also define 



POO _ „ _ 

/ 6A{u,s'H)f{uj - A) +6F{u,s''h) du 
Jo 

/■oo _ 

/ A{u;,s'H)fi^-^)du; 
Jo 



5rr"(J^) = ^ : (72) 



which gives the full fractional subleading correction to the relation between the two spectra. 
To proceed further we must introduce a model for the shape functions. 



B. Shape function models 

The shape functions are nonperturbative functions which cannot at present be calculated 
from first principles. We do, however, know several moments of these functions and we 
can use this information to constrain possible models of the shape functions. 

The leading order shape function is modeled with /mod(i^) defined in (jUUI) . We will use 
three models of the subleading shape functions. The first was introduced in based on 
the leading order function /mod('^)- The subleading functions are defined as 



^lmodl(^) — ofn 



m: 



2 J mod 



G'2modl(t^) — ~3^'^mod(^) 



-f^2modl(t^) — |/mod(^) 

mi 

imodl(^) = ^\ V modM (73) 



m; 



to reproduce the leading terms in Eq. (jS 
The second model was introduced in 
in terms of a single function 



in which the subleading functions are defined 



s{u, h) = _^e-^(^+"/^) (6(1 + Lu/A) - l) ^(1 + Lu/A). (74) 
A^ ^ 
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The dimensionless free parameter b is constrained to be 0{1) by the requirement that the 
nth moments of the functions scale hke A""*"^. We will take 6 = 1 in our plots; larger values 
of b reduce the effects of the subleading shape functions. We have 

A2 

himod2{uj,b) = -^s{uj,b) 
mi 

G'2mod2(t^,&) = -^—^S{uj,b) 

A2 

mi 

tmod2(t^,&) = --^^^—^^^5(0;, 6). (75) 

mi 

Note that in the first model the subleading shape functions vanish at u; = A, while in the 
second they are finite but nonzero. 

In our third modeF, we use a model for the subleading shape functions that has an 
additional sign flip in the region of integration. We take 

A2 
mi 



6m^ 



A2 

-f^2mod3(^) = 2/2(^) 



ml 



mi 



where 

^ ^ TT^A ^ ^ \ AJ duj \\ A/ / 
We plot the function hi{uj) in each of these models in Fig. |3] 

Although the models have very different behaviour, we can verify that they are all rea- 
sonable by calculating the first few moments Ain = S'^oo^^^'^i'^)'^^ each model and 
showing that they are of order (Aqcd/''^6)"'''^- For m^ = 4.8 GeV, for model 1 we find 
\Mi\ = (0.35GeV/mb)MA<2| = 0, {M^l = (0.35 GeV/mb)^ and {M^l = (0.29 GeV/mfe)^. 
For model 2 we find \Mi\ = (0.35 GeV/m;,)MA42| = (0.49 GeV/mb)3, |A^3| = 
(0.70 GeV/mb)^ and \A4a\ = (0.90 GeV/mft)^, while for model 3 the corresponding mo- 
ments are \Mi\ = (0.35 GeV /m^f, \M2\ = (0.54 GeV/mfe)3, IMs] = (0.54 GeV/m;,)^ and 



We thank C. Bauer for discussions of this model. 
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FIG. 4: Three models of hi(uj): model 1 (solid curve), model 2 (dashed curve) and model 3 
(dot-dashed curve). 

= (0.58 GeV/mfo)^. Thus, the first few moments of each model scale like typical 
hadronic scales to the appropriate power. Similar results are obtained for the moments of 
G2{u!) and hi{u). 




C. Numerical results 



Both the Wilson coefficients and models for the shape functions depend on the b quark 
mass rrib- While in our formulas we are implicitly using the pole mass, it is well-known that 
this leads to badly behaved perturbative series, and so we expect that radiative corrections 
to these results will be minimized if a sensible short- distance mass is used instead. The 
MS mass rfihirhh) is well-defined, but does not lead to small perturbative corrections in B 



decays j23. 



24| . The "threshold" masses, including the IS mass, PS mass and kinetic mass. 



are preferable in this context. At two-loops, a pole mass of rrih = 4.8 GeV corresponds to 
a kinetic mass m^™(l GeV) of about 4.6 GeV, PS and IS masses of about 4.7 GeV and an 
MS mass rhb{rhf)) of about 4.3 GeV. Thus, to give an estimate of the rrib dependence of our 
results, we plot them for rub = 4.8 GeV and = 4.5 GeV. 

In Fig. we plot the hadronic invariant mass spectrum using the three models of the 
previous section for the subleading corrections. These corrections are clearly large and model 
dependent over much of the spectrum. However, the integrated rate is much less sensitive 
to the subleading corrections. The functions 5Fs^(s^) and 5F^'^^(s^) defined in Eqs. f|7T|l 
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FIG. 5: Model calculations of the hadronic invariant mass spectrum dT/dsn, for rrih = 4.8 GeV 
(a) and rrih = 4.5 GeV (b). The dotted curve is the leading order result; the other curves are the 
results in the models discussed in the text. The curves are denoted as in Fig.|^ The right vertical 
line denotes the kinematic limit sh = it^'d', the left line denotes the BABAR cut sh = (1.55 GeV)^. 

and (j72p are plotted in Fig. El for the three models presented in the previous section. 

From these figures it is clear that, at least for the particular models we have chosen, the 
subleading shape functions do not contribute a large uncertainty in the extraction of \Vub\, 
and that the dominant subleading effects are from the kinematic terms. This should not be 
surprising: since there are no 0(AQCD/"^b) corrections to the total semileptonic decay rate 
Q], the subleading corrections must vanish when integrated over the full spectrum. Since 
the experimental cuts include a large fraction of the rate, the contribution to the integrated 
rate from the subleading corrections is correspondingly suppressed. This is evident from the 
plots in Fig. where the fractional correction tends to zero as the cut is increased. 

It is useful to compare these results with analogous results for the lepton energy spectrum 
in semileptonic B decays, given in ^3]- In this case, only ~ 10% of the rate is included, and 
the subleading corrections are substantial. The analogous relation to Eq. (|7(Hl is 

1/2 



where 



cfr|2\ 1/2 



rmB/2 dV jfp 



(78) 



/o""^ {niB - 2E^ - uj)f{uj - A) du 
In Fig. 13 we plot STei{E£) for nib = 4.8 GeV and rrib = 4.5 GeV in the three models used in 
this paper. It is clear from the figures that for lepton cuts near the kinematic limit Ei = 2.3 
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FIG. 6: Model calculations of the fractional corrections (JF^^' and (JF^^ to the cut width, as defined 
in Eqs. and ((7^ . for ra^ = 4.8 GeV (a,b) and rrib = 4.5 GeV (c,d). The three curves refer to 
the three different models in Fig. ^ f^F^^' includes all the subleading corrections, including those 
proportional to the leading order shape function, while STgfj^ only includes the corrections from 
subleading shape functions. The right vertical line denotes the kinematic limit sh = fni)', the left 
line denotes the BABAR cut sjy = (1.55 GeV) 2. 



GeV, the uncertainty in \Vub\ from higher order shape functions is much greater for the 
lepton energy spectrum than from the hadronic invariant mass spectrum. 



V. CONCLUSIONS 

We have calculated the hadronic invariant mass spectrum for B — > Xuii^e in terms of shape 
functions to subleading order. Introducing some simple models for the shape functions we 
have studied the dT/dsn spectrum numerically. 

Since we know little about the form of the subleading shape functions, it is difficult 
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FIG. 7: Model calculations of the fractional corrections STei, for the semileptonic B decay 
width with a charged lepton energy cut, as defined in Eq. H79() . The three curves refer to the three 
different models in Fig. ^ The solid line denotes the kinematic upper limit from B decay. 

to estimate the corresponding theoretical uncertainty in \Vuh\- However, using the spread 
of models as a guide, we can conclude that the largest subleading effects are proportional 
to the leading order shape function, and so, given a determination of the shape function 
from B Xg'-y decay, do not increase the theoretical uncertainty. Assuming our spread of 
models provides a reasonable measure of the theoretical uncertainty, we can conclude that 
the theoretical uncertainty in \ Vub\ due to higher order shape functions is at the few percent 
level. This is substantially less than the corresponding uncertainty in the integrated lepton 
energy spectrum with the current experimental cuts. This is also much less than the other 
sources of experimental and theoretical error in the current measurements of the integrated 
hadronic energy spectrum. 
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